Signatures of spin-charge separation in scanning probe microscopy 
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We analyze the effect of an auxihary scatterer, such as the potential of a scanning tip, on the 
conductance of an interacting one-dimensional electron system. We find that the differential con- 
ductance for tunneling into the end of a semi-infinite quantum wire reflects the separation of the 
elementary excitations into spin and charge modes. The separation is revealed as a specific pattern 
in the dependence of the conductance on bias and on the position of the scatterer. 



Interaction between electrons has a profound effect on 
the properties of one-dimensional electron systems. Ba- 
sic notions of Fermi liquid theory are no longer applica- 
ble, as is manifested, e.g., in the absence of well-defined 
fermionic quasiparticles related to the non-interacting 
state. Instead, the low-energy physics of these systems is 
described by the concept of a Luttinger liquid Ij. 

One consequence of the inter-electron interaction is the 
separation of elementary excitations into two branches 
with different velocities, corresponding to spin and charge 
excitations. This spin-charge separation is apparent al- 
ready within the random-phase approximation treatment 
of the interaction. The relation between the bare elec- 
tron degrees of freedom and the elementary excitations 
of a Luttinger liquid, manifested in the suppression of 
the electron density of states at low energies 2] , requires 
more sophisticated techniques such as bosonization Q . 

The suppression of the density of states was demon- 
strated in several tunneling experiments |^ Q - In con- 
trast, the separation into spin and charge modes has 
proved harder to observe. The simplest idea, to mea- 
sure the charge- and spin-density response functions, is 
problematic because coupling to the spin mode is diffi- 
cult. Another approach is provided by electron tunneling 
experiments, either into a nearby two-dimensional elec- 
tron gas (2DEG) [p, , or into a second one-dimensional 
system 0, 0, 0, linl |. The latter experiment was recently 
carried out by Auslaender et al. in quantum wires grown 
in a GaAs heterostructure In these experiments, 

the differential conductance was measured as a function 
of voltage and applied magnetic field. Patterns of os- 
cillations in the measured conductance, arising from fi- 
nite size effects, showed modulations which provide evi- 
dence for spin-charge separation in the two wires . An- 
other possibility is to use scanning tunneling microscopy 
(STM) to measure the local density of states, close to the 
edge of the s yste m , in a finite system , or next to 
an impurity |l3|. A difficulty in STM experiments is the 
need to maintain a small distance between the wire and 
the tip, necessary for electron tunneling. This is possible 
in carbon nanotubes; in recent experiments [l^ . however, 
interaction between the electrons was suppressed by the 
metallic substrate, and the results were adequately inter- 
preted in terms of non- interacting electrons. 



A scanning tip may also be placed further away from 
the wire, thus prohibiting tunneling, but still creating 
a local potential, which scatters electrons. The conduc- 
tance of the system may then be measured as a func- 
tion of tip position. This scanning probe microscopy 
(SPM) technique was successfully used for imaging the 
flow of electrons through a quantum point contact in 
a 2DEG Similar techniques were also used to 

image other systems, including confined electrons in a 
2DEG 16], multichannel one-dimensional wires 0], and 
quantum dots in nanotubes [l^ . 

In this Letter, we analyze the use of SPM in one dimen- 
sion. We consider electrons tunneling into a semi-infinite 
one-dimensional system through a fixed tunnel junction 
at its end. The scanning tip is used to create a weak po- 
tential scatterer along the wire at a distance xq away from 
its end. The tunneling conductance is then calculated as 
a function of tip position xq. Note the difference from 
the STM setup (Ref. |ilElE3) 

in which the tunneling 

occurs at the tip position. 

We find that the tunneling conductance exhibits oscil- 
lations as a function of the applied voltage V and the po- 
sition xq of the scatterer. The nature of these oscillations 
depends on the strength of the electron-electron interac- 
tion. At moderate or weak interaction, spin-charge sepa- 
ration is readily identifiable as a beating pattern, which 
allows one to extract the velocities of spin and charge 
modes. Notably, this information can be extracted from 
the voltage dependence alone, which may be useful in 
some experimental setups with a fixed position of the 
scatterer. At a strong interaction, the dominant contri- 
bution oscillates only as a function of xq, with wavevector 
2kF (here kp is the electron Fermi wavevector). 

To analyze the effect of an auxiliary scatterer on the 
current / of electrons which tunnel into the end of a one- 
dimensional wire, we express the current, 

/ cx J de[np{e) - nF{e + eV)]iyicad{e + eV) iy{e), (1) 

in terms of the densities of states i^ioad(e) and t^(e) of 
the lead and semi-infinite wire, respectively. Assuming 
t'icad(e) is essentially energy- independent, the differential 
conductance dl/dV provides a direct measure of the tun- 
neling density of states at the edge of the wire. The latter 



2 



is given by 
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where Gf'{x^x\e) is the retarded Green function, 

Gf{x,x',e)^~i die"*({V'.(x,i),^l(x',0)}). (3) 
Jo 

Here, we have taken the edge of the wire at a; = 0, and 
assumed Dirichlet boundary condition, ip{0) — 0. We 
calculate z/(e) in the presence of a weak scatterer at xq, 
as a perturbation to the known result for an ideal semi- 
infinite Luttinger liquid 

Following the standard Luttinger liquid derivation, we 
bosonize the fermionic operators, 
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Here, are the bosonic fields corresponding 

to the right (left) movers, Fg is the Klein factor, and a 
is a short-distance cutoff. The fields 3>±,s are decoupled 
into charge (p) and spin (cr) modes, as well as into their 
symmetric and anti-symmetric combinations "i? and tp, 

$±,,(a;) = ^{^p(a;) ± fp{x) + s[d„{x) ± ip„{x)]}. (5) 

At the barrier, the fields obey the boundary condition 
dxi}p{Q) = dx'&a{0) = 0, corresponding to the vanishing 
of the current at the origin. The beauty of the Luttinger 
liquid formalism is that in terms of the bosonic variables, 
the Hamiltonian becomes quadratic, 
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where Vp [va] is the velocity of the charge (spin) mode, 
and Qp, Qa are the Luttinger liquid parameters. Here, 
we neglect the interaction in the spin channel. This sim- 
plification is possible if the 2kF component of the inter- 
electron interaction potential is weak, or the energy e 
is small, as the interaction term is irrelevant under the 
renormalization fiow. In addition, in both cases = ^■ 
To calculate the tunneling density of states v{e) in per- 
turbation theory, Eqs. @-® are expressed in terms of 
time-ordered (and anti-time ordered) Green functions. 
Using the bosonization formula, Eq. the tunneling 
density of states is then expressed as 
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[Qs{t)-Qs{t)\. (7) 



Here, Qs{t) is the time-ordered Green function, 

Qs{t) - (8) 

/jigj[i?p(0,t)+siJ„(0,t)]/V2g-i[i5p(0,0)+sf„(0,0)]/V2\ 

where T is the time-ordering operator. A similar expres- 
sion holds for the anti-time ordered function Q . 

In absence of the scatterer, the average in Eq. © is cal- 
culated using the relation (e^e^)o = (^{ab+a^ I2^B^ I2)a 
(valid for A and B linear in the fields). Next, the corre- 
lation functions are calculated by expressing the fields in 
terms of the eigenmodes of Hq, 
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I'^^j2ii/p. We then have (at T = 0) 
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leading to the known power law dependence of the edge 
tunnehng density of states [l9|. 
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Here, we have set gp = g and g„ — \. 

The leading correction to v{€) due to an obstacle placed 
at X = zq, is caused by the electron backscattering from 
it. The backscattering term in the Hamiltonian, which is 
the one relevant to the calculation, is given by 
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(13) 

where U is the 2kp component of the potential created 
by the scatterer. We assume U is small, and treat it per- 
turbatively. (The opposite limit of a strong scatterer cor- 
responds to the problem of resonant tunneling through 
a double barrier, previously studied in detail j20l |: this 
limit is not suitable for SPM measurements.) 

To linear order in the scatterer strength U , one ex- 
pands the scattering matrix implicit in the definition of 
the time-ordered Green function (both in the numerator 
and in the denominator). The correction to the Green 
function ||SJ) is then given by 



dt' /2'e*[''f(°'*'+'*''"(°'*)l/^e~*[''''^°'°^+''''"^°'°^l/'^ [Hi{t') - {Hi{t'))^ 



(14) 
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Using (e^e^e'^)o = Q{AB+AC+BC+iA^+B^+c^)/2)o ^ ^j^g ^j.^^ order correction to the Green function is 
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a + ia;oSgn(< — t') + ivy\t — t'\ 
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The calculation of the tunneHng density of states, Eq. ||7J), 
thus involves integrations over t and t', along contours 
which pass several branch cuts in the complex t and t' 
planes. Here, we consider the contribution arising from 
these branch cuts assuming they are well separated, so 
that the contribution of each branch cut may be ac- 
counted for independently. The conditions for the va- 
lidity of this approximation are 
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» 1, 



exo exQ 



> 1. 
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The last of these three conditions leads to a separation of 
the contributions of the charge and spin modes. Because 
of it, the results are non-perturbative in the interaction 
strength, even in the limit of weak interaction (l—g<S^ 1). 

The contribution to the density of states induced by 
the scatterer is divided in two parts. The first one comes 
from the branch cuts at t' = xq/v^ and t — t' = xq/v^, 
where the velocities and v^, each take the value of Vp 
or fo- . Setting gp = g and = 1 , this contribution is 



U/vp Vp f a \ 



1 Vp + Vg 



l/2g+g/2-l 
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where Cpp = 1, C^^ = {vglvpfl"^^ 



2v„ V^'U 2vp 



Vp + Va 



Vp + Va 



and (j)py = (1 - l/5)7r/4 + {5pp + 5yp)TT/2. Equation ^ 
contains terms oscillating with the position of the scat- 
terer, xq, with three different wavevectors: 2{kp + e/vp), 
2{kp + tjva), and 2kp + ejvp + tjvg. Physically, these 
terms arise from processes involving a right-moving ex- 
citation (either spin or charge) traveling from the edge 
of the system, then scattering at xq into a left-moving 
excitation (again, either spin or charge) and returning to 
the origin. Note the possibility of mixing the spin and 
charge excitations in such processes, leading to the oscil- 
latory term that involves both velocities. 

The second contribution arises from the branch cuts at 



i = and i! — ±xq/v^. For this contribution, we find 
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The resulting oscillation is independent of energy. Mo- 
tivated by the weakly interacting limit in the spinless 
case, this contribution may be understood as arising from 
changes to the Friedel oscillations introduced by the scat- 
terer. Note that i5zy2 decays as x^^^'^~'^^'^ , the same power 
as the decay with position of the Friedel oscillations in 
the electron density (in the absence of the scatterer) [2l| . 

To linear order in U , the tunneling density of states is 
thus a sum of a smooth term, and oscillating terms, 
Su = Sui + i5i/2, given by Eqs. (O, and ijlHIl. The 

behavior of the oscillating terms depends on the strength 
of the electron-electron interaction. When g ^ 1, corre- 
sponding to a strong repulsive interaction, Svi rapidly de- 
cays with xq and 5v2 is the dominating term; it does not 
reveal spin-charge separation. In contrast, for a weaker 
interaction the oscillations in 5vi and 5v2 decay with xq 
with similar power laws (to linear order in 1 — g, we have 
Sui (X Xq^ and (5i^2 « Xq i+f^^^'^/^). While the decay of 
Si/i with Xq is still a bit faster than that of 6i'2, the latter 
term contains a small factor r~^{l/2g — l/2) ft! {l — g)/2. 
As a result, for a broad range of xq the main contribution 
comes from Svi. In addition, in this weakly interacting 
limit, the prefactors in Eq. H17|l may be approximated by 



pa 



Cpp = 1, to obtain 



Si'i{e) cx 



1 + sin 



exo exo 



Va 



Vn 



(19) 



X sm 2kFXQ -\ 1 1- I 

Va Vn 



The different oscillations in 5vi thus form a beating pat- 
tern. Experimentally, this may provide a tool for the 
observation of spin-charge separation. 

To illustrate these ideas further, we plot 6v in Fig. ^ 
using g = 0.8 and Vp/va = 1.4. This choice of param- 
eters corresponds to a relatively weak interaction, for 
which the assumptions leading to Eq. H19|l are valid (with 
these parameters, Caa = 0.959 and Cpa = 0.964). The 
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FIG. 1: The oscillatory part of the density of states, see 
Eqs. 1171 - II18II . in arbitrary units. The parameters used are 
g = 0.8 and Vp/vcr ~ 1.4. (a) Density plot as a function of 
e and I/xq. The main feature of the moire pattern, i.e., fast 
oscillations appearing as slanted bright and dark regions with 
a slowly modulating amplitude, is captured by Eq. 1191 . The 
contribution of Eq. 11811 appears as fainter vertical lines in the 
regions where the amplitude of the main oscillations is small, 
(b) Energy dependence at kpxo = 300. 

moire pattern expressed by Eq. (|19l) is clearly visible in 
Fig. n^a). This pattern may then be used to obtain the 
spin and charge mode velocities. 

A measurement of ^{e) may probe spin-charge sepa- 
ration even in an experiment with a fixed (and maybe 
unknown) value of Xq- Indeed, Eq. ifT^ shows that the 
number of fast oscillations of Si'i(e) fitting between two 
adjacent minima of the slow amplitude modulation is 
{vp + Va)/{vp — Va)- For example, in Fig. Qlb) there are 
6 such oscillations, corresponding to Vpjv^ — 1.4. 

The SPM method considered here has connections with 
the STM technique studied in Refs. [Hill [13. The dif- 
ferences between the two methods arise from the different 
role of the scanning tip: In STM it is used to measure the 
local density of states rather than for creating a scatter- 
ing potential for the electrons. In both cases, the mea- 
sured differential conductivity exhibits oscillations as a 
function of tip position, which contain information about 



the spin and charge velocities. The difference in setup, 
however, affects the details of these oscillations and their 
dependence on the Luttinger liquid parameters. In par- 
ticular, the SPM setup does not give rise to long wave- 
length oscillations of the type predicted for STM [1^ . 

In conclusion, the theory presented in this Letter sug- 
gests a way to study spin-charge separation in a quantum 
wire by means of a scanning probe measurement. The 
clearest manifestation of the separation is in the beating 
pattern of the differential conductance as a function of 
applied bias, see Fig. ^ The SPM method may be ex- 
perimentally useful as it only requires contacting a single 
wire at its ends, unlike STM or momentum conserving 
tunneling. Finally, we note that our results, obtained 
at T = 0, also hold at finite temperature, provided the 
distance between the tunneling point and the position of 
the probe is shorter than the thermal length v„ /T . 
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